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SYSTEMS ON LIPSCHITZ DOMAINS 
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Abstract. We develop a new approach to the L p Dirichlet problem via L 2 estimates and 
reverse Holder inequalities. We apply this approach to second order elliptic systems and the 
polyharmonic equation on a bounded Lipschitz domain f2 in M n . For n > 4 and 2 — e < 
p < 2 ^_^ + £, we establish the solvability of the Dirichlet problem with boundary data in 
L p (d£l). In the case of the polyharmonic equation A^u = with I > 2, the range of p is 
sharp if 4 < n < 21 + 1. 



1. Introduction 

Let Q be a bounded Lipschitz domain in M n . Consider the system of second order 
elliptic operators (£(u)) r = —a\^DiDjU s , where Di = d/dxi and r, s = l,...,m. We 
assume that the coefficients a[j*, 1 < i,j < n, 1 < r, s < m are real constants satisfying 
the symmetry condition aJJ = and the Legendre-Hadamard ellipticity condition: 

(i-i) ^l?l 2 l^l 2 <<;^v<-l^l 2 l^l 2 , 

/i 

for some fi > and any £ G M n , ?/ G M m . The primary purpose of this paper is to study 
the Dirichlet problem 

J £(u) = in O, 
^'^ | u = f G L p (dn) on 90 and (u)* G L p (90), 

where (u)* denotes the nontangential maximal function of u, and the boundary value of 
u is taken in the sense of nontangential convergence. The following is one of main results 
of the paper. 
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Theorem 1.3. Let O be a bounded Lipschitz domain in M. n , n > 4 with connected bound- 
ary. Then there exists e > depending only on n, m, fi and fi such that, given any 
f G L p (<9fi) with 2 — e < p < 2 ^Z^ + £ ; the Dirichlet problem (1.2) has a unique solution. 
Moreover, the solution u satisfies the estimate ||(u)*|| p < C ||f|| p . 

With Theorem 1.3 at our disposal, using area integral estimates for elliptic systems 
[DKPV] and a duality argument found in [VI], we also establish the solvability of the 
regularity problem with boundary data in L p {dVt) for p in the dual range. Here Lf(<9fi) is 
the space of functions in L p (<9fi) whose first order (tangential) derivatives are also in L p . 

Theorem 1.4. Let fi be a bounded Lipschitz domain in ~R n , n > 4 with connected bound- 
ary. Then there exists S\ > depending only on n, m, fx and fi such that, given any 
f G L p (dfl) with 2( ~™ + ^ — £\ < p < 2 + si, there exists a unique u satisfying C(u) = 
in fi, u = f on dVt, and (Vu)* G L p (<9fi). Moreover, the solution u may be represented 
by a single layer potential, and we have ||(Vu)*|| p < C||Vtf|| p , where Vjf denotes the 
tangential derivatives off on <9fi. 

We remark that for Laplace's equation in Lipschitz domains, the Dirichlet and Neumann 
problems with boundary data in L p are well understood. Indeed, it has been known since 
the early 1980's that the Dirichlet problem with optimal estimate (u)* G L p is uniquely 
solvable for 2 — e < p < oo, and the Neumann problem as well as the regularity problem 
with optimal estimate (Vu)* G L p is uniquely solvable for 1 < p < 2 + e (see [Dl, D2, 
JK, VI, DK1, Kl]). For elliptic systems as well as higher order elliptic equations, the 
solvability of the L p boundary value problems were established for n > 3 and 2 — s < 
p<2 + e (see [DKV1, DKV2, FKV, F, V3, G, PV3]). This was achieved by the method 
of layer potentials. The main tool was certain Rellich-Payne-Weinberg-Necas identities. 
We mention that for Laplace's equation, the L p estimate (2 < p < oo) for the Dirichlet 
problem follows directly from the L 2 estimate and the well known maximal principle, 
by interpolation. The L p estimate (1 < p < 2) for the Neumann problem relies on the 
classical Holder estimates for solutions of second order elliptic equations of divergence 
form with bounded measurable coefficients. Similar L°° and Holder estimates, however, 
are not readily available for elliptic systems or higher-order elliptic equations on nonsmooth 
domains. 

Nevertheless, in the case n = 3, the L p boundary value problems for the optimal ranges 
of p were solved for elliptic systems [DK2, SI, S2] and higher order elliptic equations [PV1, 
PV2, PV4] (also see [MM] for systems on manifold). In particular, it was proved that 
the Dirichlet problem (1.2) is uniquely solvable for 2 — e < p < oo. The basic idea in 
[PV2] is to estimate the decay rate of the Green's function on O, using Rellich identities, 
Caccioppoli inequalities as well as the fact that the L p Dirichlet problem is solvable for 
some p < 2. If n = 3 (or 2), the decay rate obtained is fast enough to yield the Miranda- 
Agmon maximum principle H^H^oo^ < C IMIi^dfi)) from which the L p estimates follows. 
Recently in [S3, S4], we extended this approach to higher dimensions. However, instead 
of L p estimates, we were only able to establish dilation-invariant estimates in terms of the 
nontangential maximal functions in the Morrey spaces and their pre-duals, and certain 
weighted L 2 estimates with power weights. 
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In this paper we develop a new approach to the LP Dirichlet problem. The basic idea 
is as follows. First we prove a reverse Holder inequality on 90: 




for 2 < p < 2< ~™_P + £, where A r is a surface "cube" on 90, and u satisfies Cu = 
in O, (u)* e L 2 (90) and u = on A 6r . To do this, we use the inequality (u)* < 
sup^ |u| + C/i((Vu)*), where K is a compact set in O and I\ denotes the fractional 
integral of order one on 90. The restriction of p comes from the fact that (Vu)* is in L q 
only for some q > 2. Inequality (1.5) should be considered as a localization estimate in 
the sense that if the boundary values of two solutions u, v agree on A 6r , then (u — v)* 
behaves far better than either (u)* or (v)* on A r . 

Let Ma 2t , denote the Hardy-Littlewood maximal function on 90, localized to the surface 
cube A 2r . The second step in the proof of Theorem 1.3, which is motivated by a paper of 
Caffarelli and Peral [CP], is to establish a good-A type inequality by using a real variable 
argument and estimate (1.5). Indeed, for A > Ao and the unique solution of the L 2 Dirichlet 
problem Cu = in O and u = f on 90, we will show that 

(1.6) \E(AX)\ < 5\E(X)\ + \{PeA r : M A3r (|f | 2 )(P) > 7 A} |, 

where E(X) = {P e A r : M A2r (\(u)*\ 2 )(P) > A}, and A, 5 and 7 are positive constants 
with the property 5A p / 2 < 1. The desired estimate ||(u)*|| p < C ||f|| p follows readily from 

(1.6) by an integration in A. 

The proof of estimate (1.6) relies on a C alder on- Zygmund decomposition. The key 
observation is that if 

(1.7) jj—r [ \u\ 2 da<-fX and [ |(u)*| 2 da < A, 

|Ai2r| JA 12r |Ai2r| J A 12r 

then the localization estimate (1.5) implies that 

(1.8) | {PeA r : M A2r (\(uy\ 2 )(P) > AX} | < C (J^ + ^ |A r |. 

At the end of section 2, we use these techniques to formulate a theorem of independent 
interest on the L p boundedness of operators on R n . It essentially states that if a sublinear 
operator T is bounded on L 2 (R n ), and satisfies a LP localization property in the spirit of 
(1.5) for some p > 2, then it is bounded on L q (R n ) for all 2 < q < p (see Theorem 2.32). 
This theorem may be considered as a refined version of the well-known Calderon-Zygmund 
Lemma on weak type (1, 1) bounds of L 2 bounded operators. 

We point out that the main ingredients of this new approach are, (1) interior estimates 
for solutions of Cu = 0, (2) the solvability of the L 2 Dirichlet problem, (3) the solvability of 
the regularity problem with data in L\{dVt) for some q > 2. Thus our method also applies 
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to the higher order elliptic equations studied in [DKV1, V3, PV3], where the Dirichlet and 
regularity problems in L p are solved for 2 — e < p < 2 + e. As an example, we consider the 
Dirichlet problem for the polyharmonic equation: 



g on dQ, (V^ 1 ^)* G L p (<9fi), 




ON*- 1 



where £ > 2. In (1.9), a = («i,...,a n ) is a multi-index, D a = Df 1 D% 2 ■ ■ ■ D% n , and 
| a; j = «! + ••• + a n . Also N denotes the outward unit normal to O, and q N i-i = 

Theorem 1.10. Let O 6e a bounded Lipschitz domain in R n , n > 4 roif/i connected bound- 
ary. Then there exists e > depending on n, I and O swc/i i/iat if f = {f a '■ < |a| < 
£ - 2} G WM^dfi) and 5 G L p (<9fi) raf/i 2 - e < p < + e, the Dirichlet problem 

(1.9) has a unique solution. Moreover, the solution u satisfies 

(1-n) ii(v*-yru P <c I \\ g \\ p + J2 \\VtU\\ P 

[ \a\=£-2 



We remark that WA^^dft) is a Sobolev space of the so-called Whitney arrays on dfl. 
We refer the reader to [V3] or [PV3] for its definition. As we mentioned earlier, in the 
case n = 3, the L p Dirichlet problem (1.9) was solved in [PV1, PV2, PV4] for the optimal 
range 2 — e < p < 00. It was also pointed out in [PV3] that estimate (1.11) fails in general 
for p > 2 ^J 3 1 - ) and 4 < n < 2£ + 1, or p > and n > 2£ + 2. Thus the range for p in 
Theorem 1.10 is sharp in the case 4 < n < 2£ + 1. 

The paper is organized as follows. Theorem 1.3 is proved in section 2. The proof of 
Theorem 1.10, which is very similar to that of Theorem 1.3, is sketched in section 3. Finally 
we give the proof of Theorem 1.4 in section 4. 

2. The Dirichlet Problem 

Let O be a bounded Lipschitz domain in R n . Assume G <9fi and 

(2.1) O fl S(0, ro) = {(x', x n ) G R n : x n > ifj(x')} n B(0, r ), 

where 5(0, ro) denotes the ball centered at with radius ro, ip '■ R n_1 — > R is a Lipschitz 
function. For r > 0, we let 

A r = {(i'i(i'))Gi": |xi| <r, ...,|x n _i| < r}, 

D r = {(x f , x n ) G R n : |xi| < r, . . ., |x„_i| < r, ip(x') < x n < tp(x') + r}. 
Note that if < r < cr , A r C dQ and D r C O. 
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Lemma 2.3. Let < r < ctq. There exists e > depending only on n, m, \i and O, such 
that if Cu = in Q, (u)* G L 2 (dQ) and u = on A 8r , then (u)* G L p (A r ) and 

(2A) fe/, iw H" fsc (^l, l(ur|!i ) 1 ' ! ' 

where 2 < p < 2{n — l)/(n — 3) + e and \A\ denotes the surface measure of A. 
Proof. Recall that the nontangential maximal function of u is defined by 
(2.5) (u)*(P) = sup{|u(x)| : x G O and x G j(P) } for P G dfl, 

where 7(F) = {x G O : |x - P| < 2 dist(x, <90) }. Let 

.Mi(u)(P) = sup {|u(x)| : x G 7(F) and |x-P|<cr}, 

(2.6) 

A^2(u)(P) = sup {|u(x)| : x G 7(P) and \x — P\> cr). 

Then (u)* = max { A4i(u), A^u) }. We first estimate ^2(11). To do this, we use interior 
estimates to obtain 

(2-7) |u(x)| < £ / |u(y)|dj/ < / \(u)*\do-, 

r JB(x,cr) r JA 5r 

where x G 7(P), P G A r and \x — P\ > cr. It follows that for any p > 2, 



( TaS / I^W^) " < sup.M 2 (u) 
(2.8) Vl H Ar J Ar 

- [ \(u)*\da<C [ |(u)*| 2 aV 

r\ J A 5r \|^5r| J A 5r 



< 



\A 

C 



|A 5 

Next, to estimate Mi(n) on A r , we write 



1/2 



f Xn du 

(2.9) u(x',x n ) - u(x',x n ) = - —(x',s)ds. 

Let K = x n ) G O : < cr, ^(^O +cr < i„ < cr}. Using (2.9) and interior 

estimates for Vu, it is not hard to show that for P G A r , 

(2.10) M 1 (u)(P)<sup\u\ + C f 1 p^^iS da(Q), 

K J A 2r \* *°6\ 

where (Vu)^,^ denotes the nontangential maximal function of Vu with respect to the 
Lipschitz sub-domain D 2r - By the fractional integral estimate, this implies that 

(2.11) (j^ijf |Mi(u)| p ^ /P <sup|u|+Cr(^^^ \(VvL)* D2 fdc^ '\ 
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where | = - — and 1 < q < n — 1. We now choose q > 2, depending only on n, m, 
p and so that the regularity problem for Cu = in D pr for 2 < p < 3 with boundary 
data in L\ is uniquely solvable [G] . It follows that 

(2.12) f \(Vu)* D2r \ q da < C [ \V t n\ q da<C [ \Vu\ q da, 

JA 2r JdD pr JnndD pr 

where we have used the assumption that u = on A6 r . Integrating both sides of (2.12) in 
p £ (2, 3) yields that 

(2 - i3) \t\ Jj^^ w\ L^- 

By Caccioppoli's inequality as well as its well known consequence, the higher integrability 
for Vu [Gi], we obtain 



(2.14) 



1 /•,_,„. \ /y „ / 1 



\Dzr\ J Di r J \\D± r \JD 



Vu\ q dx <C — — / iVurda; 



r 



3r 

(jD^\L 5r |u| dx ) - ~ iw^\L 5r l(u) 1 da ) 



'5r| JD 5r 

Note that g > 2 implies p > 2(n - l)/(n - 3). In view of (2.11), (2.13) and (2.14), we have 
proved that for some p > 2(n — l)/(n — 3), 



(2.15) 



1 



5r 

1/2 



-^^j' I(U) *' d ° 



This, together with (2.8), gives the desired estimate (2.4). 

To prove Theorem 1.3, it will be convenient to work with "cubes" on <9fi. Let D = 
{(x',x n ) E R n : x n > ip(x')} and dD = {(x' ,ip(x')) : x' G M n_1 }. Define the map 
$ : dD -> by VO^O) = a'- We sa Y that Q C dD is a c«6e in 3D, if $(Q) is a 

cube in M n_1 with sides parallel to the coordinate planes. A cube Q in dD is said to be a 
dyadic subcube of Q' if $(Q) is a dyadic subcube of $>(Q') in R n_1 , that is, &(Q) is one of 
the cubes obtained by bisecting the sides of &(Q') a finite number of times. Also, if &(Q) 
is a cube in R n_1 , p&(Q) denote the cube which has the same center, but p times the side 
length of &(Q). For a cube Q on dD, we denote $ _1 [p$(Q)] by pQ. As an example, in 
(2.2), A pr = pA r . 

For cube Q on <9.D and a function / defined on Q, we define a localized Hardy-Littlewood 
maximal function Mq by 

(2-16) M Q (f)(P)= sup -L / l/l^. 

Q'3P IV I JQ' 
Q'cQ 

In the next lemma, we will estimate |£?(A)|, where 

(2.17) E(\) = {PeA r :M A2r (\(uy\ 2 )(P)>\}, < r < cr Q . 
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Lemma 2.18. Suppose that 2 < p < 2 ^Z^ + £ > where e > is the same as in Lemma 
2.3. There exist positive constants A, 5, 7 and Co depending only on n, m, fx, p and Q 
such that if Cu = in O, (u)* G L 2 (dQ) and u = f G L 2 (<9fi) on <9fi, then 

(2.19) \E(AX)\ < 5\E(X)\ + \{P G A r : M Aar (|f | 2 )(P) > 7 A}| 
for all X > Xo, where 

(2.20) A = T ^/ \(uy\ 2 da. 

I^2r| JA 2r 

Most importantly, the constants A, 5 satisfy the condition 5A p / 2 < 1. 

Proof. We use a real variable argument which is motivated by the method of approximation 
in [CP]. We begin by fixing p so that 2 < p < 2(n — l)/(n — 3) + e. Let 5 G (0, 1) be a small 
constant to be determined. By the L 2 boundedness of the Hardy-Littlewood maximal 
function, 

(2.21) \E(X)\ < ° {n ' ^ l|oo) / \(uy\ 2 da. 

X JA 2r 

Thus, if A > A where A is given in (2.20) with a large C , we have |.E(A)| < 5|A r |. Let 
A = l/(25 2 / p ). Note that SA p / 2 = 1/2 P ' 2 < 1/2. 

Since E(X) is open relative to A r , there exists a collection of disjoint dyadic subcubes 
{Qk} of A r such that E(X) = [j k Qk- Clearly we may assume that each Q k is maximal in 
the sense that Q k is not contained in E{X), where Q k denotes the dyadic "parent" of Q k: 
i.e., Q k is one of 2 n ~ 1 cubes obtained by subdividing Q k once. Since l-E'(A)] < 5|A r |, we 
may also assume that 32Q k C A 2r by taking 5 sufficient small. 

We claim that it is possible to choose positive constants 5, 7 and C so that if {P G 
Q k : M A2r (|f| 2 )(P)< 7 A}^0, then 

(2.22) \E(AX) n Q k \ < S\Q k \. 

Estimate (2.19) follows from (2.22) by summation. 

It remains to prove the claim. Suppose that {P G Q k : MA 2r (|f \ 2 )(P) < 7A} ^ 0. 
Since Q k is maximal, by a simple geometric observation, we obtain 

(2.23) M A2? .(|(u)*| 2 )(P)<max(M 2Qfe (|(u)*| 2 )(P),C 1 A) 

for any P G Q k , where C\ depends only on n and HVV'Hoo- We may assume that A > C\ 
by making S small. It follows that 

(2.24) \Q k nE(AX)\<\{PeQ k : M 2Qfe (|(u)*| 2 )(P) > AX} |. 
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Also note that, if Q k C Q C A 2r , then 

(2.25) ^/ Q |f|^< 7 A »d ^J Q \W\>*<X. 

Let v = Vfc be the unique solution of the L 2 Dirichlet problem on O with boundary 
data fxi6Q fc - It follows from (2.24) and Lemma 2.3 that 



\Q k nE(AX)\ < | \peQ k : M 2Qfc (|(u- v)*| 2 )(P) > ^} 
+ |{peQ.:M 2Qfe (|(vr| 2 )(P)>^J| 

<n^rJ l(u-v)1^ + ^/ l(v)1 2 dtT 



(w/o Q J (U - Vn2 ^) P/2 + w/ 6Q J f|2d i 



where p < p < 2(n — l)/(n — 3) + e and we also used the L 2 estimate for v in the last 
inequality. This, together with (2.25), the L 2 estimate and the choice A = l/(25) 2 / p , gives 

(2.26) \Q k n E(AX)\ < \Q k \ j-^ + ^| = <5|Q fc | {c^ 2 ^" 1 + 2C275I" 1 } , 

where C 2 depends only on n, m, p and Q. 

Finally, since p > p, we may choose 5 > so small that C 2 2 p / 2 5p~ 1 < 1/2. With 5 fixed, 

we choose 7 > so small that 2C 2 ^5p~ 1 < 1/2. We then obtain \Q k n -E(A\)| < 5|Q fc |- 
The proof is complete. 

Estimate (2.19) is a good-A type inequality. Its relation with the L p norm estimate is 
well known. 

Lemma 2.27. Suppose that 2 < p < + £ where e > is the same as in Lemma 

2.3. If f G L p (dQ) and u is the unique solution of the L 2 Dirichlet problem with boundary 
data f, then (u)* e L p (A r ) and 



(2.28) 



1 



I A, 



i/p 

\(uy\ p da ' 



r\ J A r 

£Cl ^bl l(u),|2<i<r ) I/2+c (^bl |f|p ^ 1/P 



Proof. We begin by multiplying both sides of (2.19) by A 2 1 and integrating the resulting 
inequality in A on the interval (Ao, A). This gives 

(2.29) — — / A^" 1 |E(A)|rfA < 5 / \*- 1 \E(\)\d\ + C \f\ p da. 
A p l 2 Jax q Jxo JA 2r 
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It follows that 

\ pA pAXy p 

6} A5 _1 |£(A)|dA < C / AS _1 |£(A)|dA + C / \f\ p do 
J Jo Jo Ja 2t 

<CA||A 2r |+C / \f\ p da. 

J A 2r 

Since 5A p / 2 < 1, estimate (2.28) follows easily from the above inequality and (2.20) by 
letting A — > oo. 

We are now in a position to give the proof of Theorem 1.3. 

Proof of Theorem 1.3. The case 2 — £<p<2is already known [G]. For p > 2, since 
L p (dQ) C L 2 (<9fi) for p > 2, the uniqueness follows from the uniqueness for the Dirichlet 
problem with L 2 data. The existence as well as the estimate ||(u)*|| p < C||f|| p follows 
easily from Lemma 2.27 by covering <9fi with a finite number of coordinate patches. 

Remark 2.30. Theorem 1.3 also holds for the exterior domain fi_ = R n \ O. That is, 
given any f G L p (dQ) with 2 < p < 2(n — l)/(n — 3) + s, there exists a unique u on fi_ 
satisfying £u = in u = f on 90, (u)* G L p (90), and |u(x)| = 0(|x| 2_n ) as \x\ — > oo. 
Moreover, the solution u satisfies ||(u)*|| p < C ||f || p , where (u)* denotes the nontangential 
maximal function of u with respect to 

Remark 2.31. In the case n = 3, the techniques in this section yield the L p solvability 
of the Dirichlet problem for 2 < p < oo. 

We end this section by formulating a theorem, which is of independent interest, on the 
L p boundedness of operators on M. n . Its proof, which is given in [S5], may be carried out 
by slightly modifying the techniques used in the proof of Lemma 2.18. We point out that 
any operator with a standard Calderon-Zygmund kernel satisfies condition (2.33) for all 
p > 2. 

Theorem 2.32. Let T be a bounded sublinear operator on L 2 (M, n ). Suppose that for some 
p > 2, T satisfies the following L p localization property: 



(2.33) 




for any cube Q C MJ 1 and any L 2 function f with supp(f) C MJ 1 \ 3Q. Then T is bounded 
on L g (R n ) for any 2 < q < p. 




3. The Polyharmonic Equation 
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In this section we give the proof of Theorem 1.10, using the same line of argument as 
in the last section. We will need the following estimates for suitable solutions of A £ u = 
on O: 



(3.1) U^- L ur\\ 2 <C\\V"- L uU 



(3.2) IKvVpL < CllVtV*" 1 ^ 



for some q > 2 depending only on n, £, and the Lipschitz character of O. Both estimates 
were established in [V3] for any integer £ > 2. In the case of the biharmonic equation 
{£ = 2), estimate (3.1) was obtained earlier in [DKV1]. We also need a Caccioppoli's 
inequality, 

(3.3) / \V £ u\ 2 dx< . — s 9 9 / iV^updx, 

./£>„ {p-s) 2 r 2 J Dpr 

where < s < p < 1, it satisfies A £ it = in O, (V £ w)* G L 2 (A 2r ), and D Q -u = on A 3r 
for all < \a\ < £ — 1. Inequality (3.3) may be proved by using integration by parts and 
Poincare's inequalities for functions which vanish on part of the boundary. From (3.3) one 
may deduce the following boundary reverse Holder inequality for some q > 2 by a standard 
argument, 

1 ! ,_, ... _/ 1 ' ^ 1/2 



<3 - 4) \mk^ ,ix ) 

Lemma 3.5. There exists e > depending only on n, £ and O such that if A i u = in 
O, (V^ _1 ix)* G L 2 (<9fi) and D a u = /or all < \a\ < £ — 1 on t/ie surface cube Ag r , ^en 
(V* _1 u)* G L p (A r ) and 

(3 - 6) (^l |(v '" 1 " ) * |p *) 1/ ^ c (^Ij^" 1 ")* 12 *) 172 - 

where 2 < p < 2 ^J 3 1 - > + e. 

Proof. Using the well known interior estimates, one may prove that for P G A r , 

(v^- 1 «)*(p) 



n-2 



d<r(Q). 



With estimates (3.2), (3.4) and (3.3) at our disposal, the rest of the proof is similar to that 
of Lemma 2.1. We omit the details. 

As in (2.17), we let 

(3.7) F(A) = {PGA r : M A2r .(|(V £ - 1 ix)*| 2 )(P) > A} for0<r<cr . 
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Lemma 3.8. Suppose that 2 < p < 2 ^Z^ + e where e > is the same as in Lemma 3.5. 
There exist positive constants A, 5, 7 and Co depending only on n, £, p and O, such that 
if f = {f a : < \a\ < £ — 2} e WA 2 _ 1 (dQ), g e L 2 (<9fi) and u is the unique solution of 
the L 2 Dirichlet problem (1.9) with boundary data f , g, then 

(3.9) \F(A\)\<5\F(\)\ + \{peA r : M A2r {\g\ 2 + ]T | V t f a \ 2 }{P) > 7 a| | 
for all X > Xq, where 



(3.10) X = -^-l KV^uYfda. 



\^2r\ Ja 2j . ^ 

Most importantly, the constants A, S satisfy the condition 5A p / 2 < 1. 



Proof. The proof is similar to that of Lemma 2.18. However, in the place of v, we have to 
be more careful with the choice of polyharmonic function v = v k in O for each Q k , since 
WA 2 _ l (dQ) is a Sobolev space. Let (p = (pk is a smooth cut-off function on R n such that 
V? = 1 on 16Q k , <p = on dfi \ 17Q k , \D a p\ < C/r^ for < \a\ < i - 1, where r = r k is 
the diameter of Qfc. Let h be a polynomial of degree £ — 2 to be determined. We choose 
v = Vk to be the solution of L 2 Dirichlet problem (1.9) with boundary data 



(3.11) D«v = D«((u-h)p) = 3\(a-3V Up - DPk)Da ~ Pip 

at — 1 

for < \a\ < £ — 2 and Q N eJi = g<p on dQ. 

Write u = v + w + h. Then V^ _1, u = V i ~ 1 v + V i ~ 1 w, since h is a polynomial of degree 
£ — 2. It is easy to see that w is the solution of the 1? Dirichlet problem (1.9) with boundary 
values vanishing on 16Qk- In fact, D a w = D a ((u — h)(l — </?)) for < \a\ < £ — 2 and 

g N e^i = #(1 — <p) on 90. Thus we may apply Lemma 3.5 to w. 
To control the L 2 norm of (V i ~ 1 v)* : we note that for \a\ = £ — 2 
||V t V a ((w-/i)v?)||2 



(3.12) 



< C E {ll(V,(/ /3 -^/i))^-^ll2 + ||(/ /3 -^)V^-^|| 2 | 



0</3<a 

Finally, we let /i(x) = J2\ a \<e-2 iS: xa ■> where the constants B a are defined inductively 

by 

B a = , / ( P) f/a( /- :> ) for In I = / - 2 

'170* 

( 3 - 13 ) o 1 f V- 5, 

25 



17 ^-r / f a (P)da(P) for |a| 
|17Qfc| Ji 7 Q fe I , ^ (a-/3)! J 



|a|<^-2 

a>/3 
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for < | a | < £ — 2. We remark that B a is defined in such a way that 



/ {fp - D p h) da = for all < \(3\ < £ - 2. 

J17Qi. 



(3.14) 

'l7Q k 

From this, by using Poincare's inequality repeatedly, we obtain 

\\h ~ D P h\\ L , {17Qk) <Cr D a h\\ LHl7Qk) 

\<*\=\P\+i 



(3.15) <<7r^l- 2 £ ll/a - ^Ik-dTO.) 

|a|=£-2 

< c^-i^"- 1 £ iivj a ii L2(17Qfe) . 

|a|=^-2 

By combining estimates (3.12) with (3.15), we get 

(3.16) ||( V ^)*||2<c{|b|| i2(170fe) + l|V*/a|U-(17Q fc )j- 

^ \a\=£-2 > 

With these observations, the argument in Lemma 2.18 goes through with minor changes. 
We leave the details to the reader. 

The following lemma follows from Lemma 3.4 by integration. 

Lemma 3.17. Suppose that 2 < p < 2 ^_P + s, where e > is the same as in Lemma 
3.5. If f G WA^_ t (dQ), g G L p (dQ), and u is the polyharmonic function on O satisfying 

(V' _1 u)* eL 2 (dtt) andD a u = f a forO< \a\ <£-2, =9 ondVt, then (V^u)* G 

L p (A r ) and 

i/p < i r ^ V2 



(3.18) 



As before, Theorem 1.10 is a consequence of the L 2 solvability and Lemma 3.17. 
4. The Regularity Problem for Elliptic Systems 

This section is devoted to the proof of Theorem 1.4. We follow an approach found in 
[VI] , where it was used by Verchota to establish the solvability of the L p regularity problem 
for Laplace's equation. The basic idea is to prove 

(4-1) ||_^|| p <C||V t u|U 

for star shaped Lipschitz domains, using a duality argument which involves conjugate func- 
tions and area integral estimates. With (4.1) one may show that the single layer potential, 
which maps L p (dfl) to L\(dVl), is invertible. We should mention that in (4.1), ^ denotes 
the conormal derivative defined by (f^) r = a ij^~^i- 
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Lemma 4.2. Let Q be a starshaped Lipschitz domain in M. n , n > 3. Suppose that for some 
p > 1, the L p Dirichlet problem (1.2) is uniquely solvable for any f G L p (dQ). If Cu = 
in £1, (Vu)* G L p (<9fi) and Vu has nontangential limit on dVL, then ||§^|| P ' < C ||Vtu|| p > 
where p' = p/(p — 1) . 

Proof. We may assume that is starshaped with respect to the origin. Let g be a Lipschitz 
continuous function on <9fi and v be the solution of £v = in O satisfying (Vv)* G L 2 (<9fi) 
and v = g on <9fi. Define 

C 1 dr 

(4.3) H(x) = / w(rx) — for x G O, 

where w(x) = v(x) — v(0). It is easy to verify that £H = and Xi = w r . As in the 
case of harmonic functions, we also have 

(4-4) ll(VH)l P <C||(w)%<C||g||^ 

The first inequality in (4.4), which holds for any p > 0, follows from the area integral 
estimates for elliptic systems [DKPV], and the second by the solvability of the L p Dirichlet 
problem. We refer the reader to [V2, p870-871] for details in the case of Laplace's equation. 

We now use integration by parts to obtain 

f du f 9v , f dw , 

/ —-gda= u-—da= u-—da 

Jon Jdii diy J m dv 



(4.5) 



. f dH f du r rs dH s 

= (2 - n) I u • — da- j^ a zj ~^ — x e da, 



where = rii-^ — n£-^r- is a tangential derivative. We remark that the integration by 
parts used in (4.5) can be justified by an approximation argument, using estimate (4.4) 
and the assumption that Vu has nontangential limits on dtt and (Vu)* G L p . It follows 
from (4.5) and Holder's inequality that 



(4.6) 



-^-■gda 



<C||V t u|U|g| 



pi 



where (4.4) and Poincare's inequality are also used. The desired estimate then follows by 
duality. 

We now give the proof of Theorem 1.4. 

Proof of Theorem 1.4. Let T(x) denote the matrix fundamental solution for the oper- 
ator C on M n with pole at the origin. Consider the single layer potential 

(4.7) S{g){x) = [ T{x-y)g{y)da. 

JdCl 
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For the existence as well as the estimate ||(Vu)*|| p < C|| V t u|| p , it suffices to show that 
there exists e > such that S : L p (dfl) — > L\(dQ) is invertible for 2 ^~ L 1 - ) — e < p < 2. To 
this end, we fix P G <90 and assume that 

(4.8) B(P ,r )nn = B(P Q ,r Q )n{(x',x n ) Gl™ : x n > tfj(x')} 
where ip : R n_1 — > R is Lipschitz continuous. We may assume that Po = 0. Let 

(4.9) r = {(i',i n )6l": < r, . . . , |x n _i| < r, ip(x') < x n < Br} 

where the constant B = B(n, ||V^||oo) > is chosen so that fi r is a starshaped Lipschitz 
domain for any r > 0. Note that, by Theorem 1.3, there exists 6 > depending only on n, 
m, [i and ||VV>||oo) such that the L q Dirichlet problem for the operator C on Vt r is uniquely 
solvable for any r > and 2 < q < q = 2 ^Z^ + <5- 

Let g G L p (dQ) and v = S(g) in R n . Suppose 2 < p' < qo- We may apply Lemma 4.2 
to v on O r for < r < cr . This gives 

(4.10) / \^\ P da<C [ \V t v\ p da<C [ \V t v\ p da + C [ |Vv| p da, 

where s G (1,2). By integrating (4.10) in s and covering dfl with a finite number of 
coordinate patches, we obtain 



(4.H) 



\^\ P da<C [ \V t v\ p da + C f |Vv| p dx 



<C [ |Vtv| p da + 7 / |(Vv)*| p rfa + C 7 sup|v| p , 



where O = {x G O : dist(x, dCl) < r }, and K is a compact set in O. To estimate v on K, 
we choose g < 2 so that the LP Dirichlet problem is uniquely solvable. It follows that 

(4-12) sup |v| < C ||(v)*|| 5 < C \\v\\ q < C{\\V t v\\ p + ||v||p}, 

where | = i — and we used the Sobolev imbedding in the last inequality. Note that 
q < 2 implies p < 2(n — l)/(n+ 1). We may assume that p' < 2(n — l)/(n — 3) + 5. Thus, 
if p < p < 2, we obtain 

(4.13) j \^\ P da<C [ \V t v\ p da + C [ |v| p da + 7 C / |g| p rfa, 
Jan 0^ ian Jan Jan 

where we have used the estimate ||(Vv)*|| p < C ||g|| p (see e.g. [VI] or [Kl]). By the same 
argument, estimate (4.13) also holds for the exterior domain R n \ O. Thus, by the jump 
relation g = — where ± indicate the limits taken from O and R n \ Q respectively, 
we obtain 

(4.14) ||g|| p < ||-^±|| p + || — || p < C {||Vv t || p + ||v|| p + 7 ||g|| p }, 
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for p < p < 2. Choose 7 > small so that 7C < 1/2. This gives 
(4-15) ||g||p<C||<S(g)|| Lf(an) . 

Estimate (4.15) implies that the operator S : L p (dQ) — > L\(d£l) is one-to-one and has a 
closed range. Note that the range is also dense, since S : L 2 (<9fi) — > L\ (<9fi) is invertible 
[G]. We conclude that S : L p (dQ) -> L?(dfi) is invertible for p < p < 2. 

Finally, to show the uniqueness, we suppose that £u = in O, u = on dfl and 
(Vu)* G L p for somep > 2(n - l)/(n + 1) -e. It follows from (2.10) that (u)* G L q , where 

- = 737- Note that if e > is small, q > 2 — e. Thus the uniqueness in Theorem 1.4 

follows from the uniqueness of the L q Dirichlet problem for q > 2 — e [G]. 
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